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{Article by Guo Bijun [6665 3880 0689], Deng Longjiang 
(6772 7893 3068], and Chen Xuanzhen [7115 6693 
5271] of the Department of Materials Science and Engi- 
neering, UEST, Chengdu 610054: “Hexagonal Ferrite 
Ku-Band Radar Absorbers”; MS received 8 Sep 91, 
revised 28 Oct 91; project supported by the Preliminary 
Research Fund of MMEI] 


[Text] Abstract 


The design concept and preparation technique of hexag- 
onal-ferrite-based Ku-band radar absorbent materials 
(RAM) is discussed. The physical properties and mea- 
sured electromagnetic parameters are presented. Finally, 
preliminary application results are described. The absor- 
bence is greater than 12 dB in the 2-cm band. If the 
composition of the coating is further adjusted, an absor- 
bence of greater than 10 dB can be achieved in the 2-3 
cm band (8-18 GHz) range. 


Key Words: Ku-band, ferrite, absorber, spontaneous 
resonance, radar absorber, coprecipitation technique 


Classification: TM277 


I. Introduction 


The use of magnetic materials to absorb microwave energy 
is based on spontaneous resonance. When the microwave 
frequency coincides with the natural resonance frequency of 
the material, maximum absorption occurs. The spontaneous 
resonance frequency is proportional to the anisotropic mag- 
netic field inside the material. Hexagonal ferrites have a 
high anisotropic field. Therefore, their natural resonance 
frequencies are high and can be used to absorb microwave 
energy. Elements such as Co and Ti are used to replace the 
Fe*** ion in BaFeO,, to produce a microwave absorber 
Ba(CoTi),Fe,>,,0,;, by ion substitution. The anisotropic 
field can vary over a wide range as a function of x. 
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Figure 1. Ba(CoTi), Fe,,..0,. Natural Resonance Fre- 
quency vs. x 











Consequently, it 1s possible to produce a wide-band micro- 
wave absorber. Figure | shows how the natural frequency of 
Ba(CoTi),Fe,>.O,, decreases as a function of x 


In addition to composition, the shape, size and orienta- 
tion of the absorbent power have significant effects on its 
absorption characteristics. When the power size 1s 
reduced to the monodomain level, the domain-wall 
resonance mechanism no longer exists. At this time, the 
natural resonance frequency of the power increases 
slightly. However, the maximum microwave energy it 1S 
capable of absorbing falls slightly as well, as shown in 
Figure 2. The powder prepared in this work 1s approxi- 
mately 7-10 um in size 
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Figure 2. Effect of Powder Size on Microwave Absor- 
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Figure 3. Effect of Magnetic Field Orientation on 
Microwave Absorbence a 





II. Experimental Results 


High-purity BaCO,, Fe,O,, CoO and 710, are used as 
raw materials. They are weighed stoichiometrically and 
then ground in a ball mill for 8 hours. The powder is 
pre-heated at 1280-1350°C, crushed and then sintered. 
After it is crushed again and sieved, it finally goes 
through an annealing process. The physical properties 
and electromagnetic parameters of different M-phase 
powders of different x content were experimentally 
determined. The results are as follows: 











1) A scanning electron micrograph of the powder, shown 
in Figure 4, clearly proves that the powder is hexagonal 
and the average particle size isa < 10 pm. 














i = ad 
~ F 
-»” 
= 
t 


Figure 4. Scanning Electron Micrograph of a 


Hexagonal Ferrite Ku-Band Radar Absorber 





2) The X-ray diffraction analysis result, shown in Figure 5, 


indicates that the material is a single M-phase material. 


3) The specific saturated magnetization moment o, is 
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Figure 5. X-Ray Diffraction Pattern of Ferrite Ku-Band 





measured using a magnetic balance and o, = 54.0 (emu/g). 


4) An HP8510B network analyzer was used to determine 
electromagnetic parameters such as relative complex 


magnetic permeability py, = p,’ - jp, and relative com- 


plex dielectric permeability €, = €,’ - je,”. The results are 
shown in Table 1. 





Table 1. Relative Complex Magnetic Permeability and Dielectric Permeability, HP8510B Network Analyzer 


























F (GHz) E,’ E,“ E,"/E;" U;’ U,* U/"/Ur 
12.40 6.955 0.076 0.011 1.796 0.323 0.100 
13.50 7.032 0.079 0.011 1.784 0.727 0.407 
14.50 7.036 0.121 0.017 1.360 0.998 0.729 
15.50 7.031 0.100 0.015 1.001 0.830 0820 
16.50 7.021 0.119 0.017 0.924 0.616 0.667 

17.50 7.021 0.119 0.017 0.954 0.520 0.546 

18.00 7.078 0.061 0.007 0.956 0.516 0.540 























III. Application 


1) The initial Ku-band RAM made have no other additives besides an epoxy binder. Table 2 shows the reflective 
characteristics of the coated specimens. 





Table 2. Reflective Characteristics of RAM Made of Ferrites and Epoxy Binder 
(absorber thickness 1.5 mm, area 180 x 180 mm/’, single layer, backed by metal plate) 














Frequency (GHz) Reflectivity (4B) Frequency (GHz) Reflectivity (4B) 
12.40 3.22 15.40 9.38 — 
12.60 -3.86 15.60 914 © 
12.80 4.62 15.00 8.62 
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Table 2. Reflective Characteristics of RAM Made of Ferrites and Epoxy Binder 
(absorber thickness 1.5 mm, area 180 x 180 mm’, single layer, backed by metal plate) (Continued) 

Frequency (GHz) Reflectivity (4B) Frequency (GHz) | Reflectivity (4B) 
13.00 6.24 16.00 -8.06 

13.20 6.98 16.20 -8.48 
13.40 -7.68 16.40 -8.74 
13.60 -8.58 16.60 -8.90 
13.80 -9.72 16.80 -9.34 
14.00 -10.50 17.00 
14.20 -10.74 17.20 
14.40 -10.88 17.40 -10.06 
14.60 -10.78 17.60 -9.86 
14.80 -10.26 17.80 | 10.54 
15.00 -9.90 18.00 -4.04 
15.20 -9.86 
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2) The reflective characteristics of specimens coated with a combination of ferrite and other absorbers are shown in 
Tables 3 and 4. 





Table 3. Reflective Characteristics of Microwave RAM Made of Ferrite and Other Absorbers 
(absorber thickness 1.29 mm, area 180 x 180 mm’, single layer, backed with metal plate) — 







































































Frequency (GHz) Reflectivity (dB) Frequency (GHz) Reflectivity (dB) 

12.40 -10.20 i640 t—t™” L 20.42 

12.60 -11.40 15.60 | 22.52 

12.80 -12.70 15.80 | 22.26 

13.00 15.60 16.00 | 24.92 

13.20 -17.02 16.20 | 24.54 

13.40 -20.60 16.40 | 2466 
13.60 24.82 16.60 | 24.42 

13.80 -26.28 16.80 T2288 

14.00 24.50 17.00 Tne? 

14.20 22.78 L 17.20 T -_ 19.52 oe 
14.40 21.02 | 17.40 17.64 

14.60 20.54 17.60 TT 610. 
14.80 -19.62 17.80 «450 
15.00 20.38 18.00 _ il 10.78 

15.20 20.86 f a 














Table 4. Reflective Characteristics of Microwave RAM Coated With Ferrites and Other Thickening Agents 
(absorber thickness 1.26 mm, area 180 x 180 mm’, single layer, backed with metal plate) 





























Frequency (GHz) Reflectivity (4B) Frequency (GHz) si Reflectivity (4B) 

12.40 12.06 15.40 a i ee 
12.60 -12.98 15.60 | -20.56 

12.80 -13.80 15.80 | 19.52 

13.00 -15.68 16.00 .206602~¢O~*~*” 
13.20 16.94 16.20 —.2042—~CO*~<CS~*S 
13.40 -17.90 16.40 [20.40 
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Table 4. Reflective Characteristics of Microwave RAM Coated With Ferrites and Other Thickening Agents 
(absorber thickness 1.26 mm, area 180 x 180 mm/?, single layer, backed with metal plate) (Continued) 
































Frequency (GHz) Reflectivity (dB) Frequency (GHz) Reflectivity (dB) 
13.60 “19.12 16.60 | 20.30 
13.80 19.72 16.80 | 18.98 
14.00 -19.20 17.00 | 18.82 
14.20 -19.34 17.20 | 17.32 
14.40 -18.56 17.40 16.06 
14.60 -18.94 17.60 T 15.32 
14.80 -18.82 17.80 T 14.34 
15.00 -19.58 18.00 | 12.52 
15.20 20.14 T 











3) Coated specimens weigh approximately 3.52 kg/m?. 


From the reflective characteristics of the two types of 
coatings shown above, we can see that: 


1) The absorbence of the specimen can reach above 12 
dB over the entire 2-cm band. This indicates that such 
absorbers can cover the entire Ku band. 


2) If the absorbent coating formula and electric structure 
can be further adjusted, it is highly probable that a 
wide-bandwidth, high-absorbence coating which covers 
the 2-3 cm (8-18 GHz) band with an absorbence of more 
than 10 dB can be made. 


3) Compared to conventional ferrite-based RAM, the 
Ku-band coating developed is thicker, but lighter in 
weight. This is especially critical for an aircraft. 


IV. Conclusions 


Using ion substitution, it is possible to control the 
microwave operating frequency by replacing some of the 
Fe*** in BaFe,,0,, The new hexagonal ferrite developed 
is an excellent Ku-band absorber. The successful devel- 
opment of such a RAM is of significant importance to 
the preparation of wide-band stealth materials. 
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[Excerpts] Abstract 


The scattering characteristics of multilayer frequency- 
selective surfaces (FSS) are analyzed by the method of 
moments in conjunction with generalized scattering 
matrix theory. Through the introduction of a generalized 
waveguide, a unified treatment of a multilayer FSS with 
an arbitrarily shaped element is realized. As examples, 
the scattering characteristics of several FSS structures 
are computed. The results are in good agreement with 
published experimental data. 


Key Words: electromagnetic wave, scattering, frequency- 
selective surface. 


I. Introduction 


It is well known that a FSS has the property of a 
bandpass or band-rejection filter. It is widely used in 
radomes, reflector antennas and beam-control systems. 
A great deal of work has been done on the analysis of 
single-layer FSSs. Many methods have been presented, 
including mode matching,':’ spectral domain” and sub- 
domain moment.® However, these methods are universal 
and have their limitations. As for the analysis of multi- 
layer FSSs, it is limited to two layers. ’"* 


In this work, the scattering characteristics of a multilayer 
FSS are analyzed using the method of moments and 
generalized scattering matrix theory. The basic idea is to 
introduce the generalized waveguide concept and con- 
sider each element on the FSS as a joint for two gener- 
alized waveguides. The intrinsic mode of a generalized 
waveguide is the periodic Floquet mode.'? The method 
of moments is used to find the generalized scattering 
matrix at the waveguide joint. Then, the scattering 
parameters of the multilayer FSS can be calculated based 
on the generalized scattering matrix. 
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II. Theoretical Analysis 


The following is an analysis of the theory using an infinite 
conductor screen with a periodic two-dimensional array of 
holes. For a two-dimensional planar bonded-plate array, 
the analysis can be derived from the results of this work on 
the basis of duality. The structure studied is shown in 
Figure |, where d, and d, are the periods along the x-axis 
and s-axis, respectively, and a is the angle between the two 
axes. [passage omitted] 
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Figure 1. Gelber of the FSS 





III. Numerical Results 


The theory discussed above is used to calculate the 
reflection and transport characteristics of a plane wave 
incident upon a single- and a double-layer FSS. Specifi- 
cally, the FSS scattering characteristics of infinite plates 
with an array of rectangular and circular holes are 
calculated in this work. In the computation, W, = M,. 
Furthermore, in order to be universally applicable, the 
intrinsic mode of the cross section of the waveguide, in 
the shape of the hole on the plate, is chosen to express 
M,. In the actual computation, Q = 10, i.e., one con- 
siders the lowest 10 modes in the waveguides. In equa- 
tion (1), p and q should be selected so that t,, < 11 k 
(where k is the wave number in free space). Figures 3 and 
4 show the reflection and transmission characteristics of 
a vertically incident TE wave on a single-layer FSS. The 
transmission frequency characteristics of a vertically 
incident TE wave on a double-layer FSS (comprised of 





rectangular arrays) are shown in Figure 5. Other results 
in existing literature are also shown in the figures. It is 
obvious that they are in excellent agreement. 
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Figure 3. Transmission Characteristics of Main Mode 
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Figure 4. Reflective Index of Main Mode From Single- 
Layer Circular Array 
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Figure 5. Power Transmission of Double-Layer Rectan- 
gular Array FSS 
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= |) 














IV. Conclusions 


The scattering characteristics of multilayer frequency- 
selective surfaces (FSS) are analyzed by the method of 
moments in conjunction with generalized scattering 
matrix theory. The theory can be applied to any multi- 
layer FSS structure with any shape of elements. The 
analysis done in this work is based on the assumption 
that the thickness of the FSS is infinitely thin. Neverthe- 
less, it can be easily extended to multilayer FSSs with a 
finite thickness. 
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[Text] Abstract 


Motion compensation is a key problem in ISAR (inverse 
synthetic aperture radar). Usually, compensation 1s 
made using an isolated point on the target. In reality, this 
isolated point may not exist. This paper presents a new 
scheme to compensate for the echo of a continuous 
target. 


Key Words: inverse synthetic aperture radar, radar 
imaging, motion compensation. 


I. Introduction 


ISAR imaging belongs to small-rotation-angle observa- 
tion imaging with a rotation angle ranging from 3°-5°. It 
is known! that targets such as airplanes, missiles and 
ships are considered as a series of discrete scattering 
points (perhaps very concentrated) as far as a high- 
resolution radar is concerned. Within an azimuth angle 
of less than 30°, these scattering points can be considered 
approximately fixed and the scattering strength is a 
constant. Hence, ISAR imaging is used to plot the 
distribution of these scattering points. If the target is on 
a rotating platform, this is not difficult to accomplish. 
Therefore, that case will not be discussed. 


A real moving target, as far as a radar is concerned, has 
a translational component in addition to rotation that 
contributes to imaging. Hence, it 1s necessary to select a 
reference point on the target prior to ISAR imaging in 
order to compensate for the translational component of 
the reference point. This is equivalent to placing the 
reference point on the rotation axis. The imaging process 
is carried out based on compensated data. 


ISAR imaging can handle a variety of signais. In this 
paper, the discussion is limited to narrow pulses. The 
results can be extended to other types of signals. 


II. Concept and Purpose of Compensation 


In ISAR imaging, there is relative motion between target 
and radar. As shown in Figure |, using the target as the 
reference, i.e., choosing a point on the target as a stl 
reference point, the radar moves along the broken line 
(with dashes and dots) shown in the sketch. As far as 
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rotating platform imaging is concerned, the radar should 
be in circular motion along the arc shown in the sketch. 
Therefore, the purpose of motion compensation is to 
provide proper adjustment to the time delay of the echo. 
(Intensity should also be adjusted. However, if relative 
distance does not vary much, intensity variation can be 
neglected.) 





Figure 1. Relative Motion Between Radar and Target 





Assume an echo is expressed as follows: 
s(t) = f(t) exp(j@ot) (1) 


where f(t) is the complex envelope of the echo and @, is 
the angular frequency of the carrier. 


By delaying the echo by a time interval t, the delayed 
signal is: 


s(t-t) = fit-t) exp(-j@_t) exp{jWot) (2) 


Usually, radar signals are processed at the fundamental 
frequency after mixing. After multiplying the above 
equation by exp(-jWot), we obtain the delayed complex 
envelope f(t-t) exp(-jw>t). When the time delay t is 
changed, the complex envelope f(t-t) and initial phase 
factor exp(-j@ot) will also change. However, the rate of 
change is different. For an X-band radar, a delay of the 
order of a millimeter causes significant change in its 
initial phase. On the other hand, it almost has no effect 
on the complex envelope when the pulse width is several 
dozen centimeters. Based on this reason, and the fact 
that initial phase changes come around every 22, motion 
compensation can be performed in two steps. The first 
step is to carry out time delay translational correction 
(rough compensation). The accuracy is of the order of a 
centimeter. The next step is to further adjust the initial 
phase of the corrected echo (fine compensation). The 
error should be far less than x/2. Upon completion of 
these two steps, the initial phase is essentially accurate. 
However, the complex envelope delay might have an 
error of the order of a centimeter. This should not have 
any significant effect on imaging. 


Hence, if the target has an echo from an apparent 
isolated scattering point, motion compensation is not 
difficult. It can be done by simply choosing that isolated 
point as a reference point to align the envelope and to 
change the initial phase so that the phased difference 
across a cycle is zero (which is equivalent to the axis of 


rotation). However, only a few special cases have iso- 
lated scattering points. Many targets have continuous 
echoes that are a superposition of echoes from numerous 
scattering points. 


In the event that there is no isolated point, a point from 
which its delay and phase shift can be determined should 
be chosen as the reference point. In the following, we 
shall prove that it is feasible to use the “scattering center 
of gravity” as a reference point. 


III. Calculation of Time Delay and Phase Shift of 
Scattering Center of Gravity 


The longitudinal and transverse coordinates of the scat- 
tering center of gravity of a group of scattering points can 
be calculated based on their positions and echo intensi- 
ties. In reality, only echo data ts available. The problem 
is how to calculate the coordinates from echo data alone. 


In order to demonstrate how the longitudinal time de!ay 
of the scattering center of gravity can be calculated from 
echo data, let us use the case shown in Figure 2 as an 
example. In that case, the echoes do not overlap on the 
time axis. In the kth observation, if the echo intensity 
and time delay of each scattering point is o, and t, (i = 
1, 2, ...), then the longitudinal coordinate T, of its 
scattering center of gravity is 


Pe= Ditie/ Do, (3) 


Electromagnetic waves 
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Figure 2. Target Scattering Points onl Echoes 





However, it is not so simple in real life. The echoes from 
the scattering points of a continuous target do overlap. In 
this case, the kth echo can be expressed as: 


$4(t)= Di /o,v(t—-t,,) (4) 


where v(t) is the transmitted signal. The complex enve- 
lope of the echo is 


f(t) = 2ivV/o,(t —t, gexp(—j Wot; ») (5) 








where a(t) is the complex envelope of v(t). It is obvious 
that 


fat) ?=fe Ct F2(t) = DJ0, |a(t—t,,) |? 


+ Di Biv9;0,4(t 84): 


i*! 
-a*(t—2t,,)exp(—J0; 12) (6) 


where Os x=O5n—Oip=Oglt 4 —t 54) 


The first term in equation (6) is the power of each 
scattering point. The second term represents the cross 
term. In practice, it is only necessary to consider the 
overlapping part of the echoes. As mentioned earlier, t,, 
varies only slightly for neighboring echoes and has no 
effect on a(t). However, phase may vary considerably in 
a random fashion. In an approximation, 8,, can be 
treated as a symmetrically and randomly distributed 
variable with a zero mean value for each observation. 
Hence, when taking the statistical average of (6), the 
cross term is zero, i.e., 


EC fx(t) "I= Doo, lace —t,4)/? (7) 


This suggests that there is no difference in the calculation 
of echo power compared to the case where all the echoes 
do not overlap from the standpoint of statistical average. 
Nevertheless, in order to calculate the scattering center 
of gravity (as shown in equation 3), it is necessary to use 


Daitie. 


The echoes from the target can be divided into range 
rings according to different side lobes. There might be 
many scattering points inside a range ring (with a spe- 
cific t,, value). During continuous observation, due to 
phase variation, the echo also fluctuates. However, as a 
statistical average, it is identical to the situation in which 
every range ring has a single scattering point. In other 
words, when calculating the time delay of the longitu- 
dinal coordinate of the scattering center of gravity with 
equation (3) based on the echoes obtained, it also fluc- 
tuates. However, its average is the correct value. 


As a verification, we designed a moving target with 20 
scattering points. Figure 3 shows the time delay of the 
scattering center of gravity calculated from simulated 
echoes. The disturbance is apparent. After smoothing, it 
is very close to the original design curve. We also added 
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some noise to the simulated echoes. The effect of noise 
amplifies the disturbance, which can be reduced by 
smoothing. 
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Figure 3. Translational Trajectory of Scattering Center 


of Gravity 





After delay compensation of the complex envelope, the 
next step is initial phase compensation. The motion of a 
target can be divided into translational and rotational. 
The periodic phase shift of the echo from each scattering 
point is the sum of the phase shifts produced by these 
two types of motion. The translational phase shift is the 
same for every scattering point, while the rotational 
phase shift might be different from point to point. With 
respect to the rotation axis (reference point), rotational 
phase shift is zero. Therefore, when the scattering center 
of gravity is chosen as the reference point, the phase shift 
at the center of gravity (not necessarily a scattering 
point) should be calculated and used for initial phase 
compensation. 


If rotation alone is taken into consideration and a 
stationary reference point is chosen as the origin, then 
the longitudinal displacement of scattering point i, Ay 
generated as a result of rotation by an angle Ag is (see 
Figure 4): 


Ay; = R, cos (9; - Ag) - R, cos 9; = y,; (cos Ag - 1) + x, sin 
Ag (8) 


Electro- 
magnetic 


wave {x (transverse) 





y (longi- 
tudinal) 








Figure 4. Longitudinal Displacement of Scattering Point 
Caused by Target Rotation 
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In reality, the target is several hundred kilometers away 
and two pulses are separated at most by a few millisec- 
onds; the Ag generated by target movement within this 
time interval is usually less than 0.01°. If x; and y, are 10 
meters, then the first and last terms of equation (8) are of 
the order of a micron and millimeter, respectively. 
Hence, 


Ay; = x; Ag (9) 


As far as the ith echo is concerned, the phase-angle 
difference between two consecutive observations due to 
target motion is 


®, = (4n/A,) x, Ao (10) 


If the transverse coordinate of each scattering point is x, 
(i = 1, 2...), then the transverse scattering center of 


gravity is 


To= 2st i/ 319; (11) 


By multiplying (42/A,) Ag to both sides of the equation 
above, and by taking equation (10) into consideration, it 
is possible to obtain the phase-angle difference of the 
transverse scattering center of gravity after the target 
rotates by an angle of Ag. 


P= D0: 951 PIF (12) 


In reality, the distribution of scattering points is 
unknown. We need to find ®, from the waveform 
observed. If it is simple as is the case shown in Figure 2 
where the echoes do not overlap, the phase shift ®, can 
be calculated from the phase change of consecutive 
echoes. Then, ®, can be obtained from equation (12) 
based on the intensity of the echo o,. 


If the scattering center of gravity is chosen as a stationary 
reference point, then the phase-angle difference between 
observations, ®,, should be zero. Using the ®, obtained 
from equation (12) as the initial phased compensation 
should produce the desired result. 


In real life, many sub-echoes may overlap at the same 
time. In this case, the phase changes of these sub-echoes 
cannot be determined. In order to obtain ®, from (12), 
the easiest way 1s to replace o, with echo power at a 
specific time (usually in an interval of a pulse width) and 
to use the phase difference of the echo at this time 
interval between two consecutive observations as ®,. 
Does this arrangement make sense? Let us take a look at 
the case where only the echoes from two scattering points 
overlap. 


As we mentioned several times earlier, in a few neigh- 
boring observation periods, as an approximation, the 
complex envelope can be considered as stationary. This 


means that the echo vector from any scattering point al 
any instance remains unchanged. Only the phase angle 
changes. Figure 5 shows the situation involving two 
scattering points in two consecutive observations. In the 
figure, r, and r, are the amplitudes of the echoes at those 
two scattering points at that instance; R, and R, are the 
synthetic amplitudes of the two observations. Using the 
vector of r, as a standard, the phase angle varies by Ay. 
The vector angles of the two sub-echoes in the two 
observations are 6 and 6 + A@, with 6 = 6, - 6, and A@ = 
A6, - A®,, where 6,, 6, and A@,, A@, are the phase angles 
of the sub-echoes and the phase-angle differences in the 
second observation, respectively. (The variation of 8, 1s 
not shown in the figure.) If A@ and Ay are both small, 
then geometrically we can determine that 


R, = R,=R (13) 
Ay = r,A8 cos (8 - y/R (14) 


If o, and ®, in equation (12) are directly replaced by the 
echo, the contribution to equation (12) at that moment is 
R7(Ay + A6,). After substituting it into equation (14) and 
taking the geometry shown in Figure 5 into account, we 
have 


R7(Ag + AO,) = (r,7 + r,r, cos 8) AO + (r,7 + 1,7 + 21,1, Cos 


8) A, (15) 
Taking the average of the above equation over 9, we get 
E[R7(Ag + A@,)] = r,? AB, + r,7 AB, (16) 


This demonstrates that the phase-angle difference of the 








Figure 5. Variation of Echo From a Rotating Target 





transverse scattering center of gravity derived from the 
power and phase difference of the measured waveform is 
consistent, from the standpoint of statistical average, 
with that calculated based on all scattering points. There- 
fore, the phase compensation calculated from the wave- 
form must be filtered to remove disturbance from it. 


The above is the situation when two sub-echoes overlap. 
It is easy to extend to more overlapping sub-echoes. In 
the above discussion, it is also assumed that A@ — 0. 
When A®@ has a finite value, direct analysis is more 
difficult. Through computer simulation, it was found 
that there is a difference between the statistical average 
of the calculated transverse scattering center of gravity 
and the real value. However, as long as the absolute 
value of A® < 0.5 rad (which covers most cases), the 
error is less than 5 percent. This kind of error is tolerable 
in ISAR imaging. 
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With a repetition frequency of the order of a few 
milliseconds, the phase-angle differences of various scat- 
tering points caused by rotation are small. However, 
translational motion will cause a larger but equal phase 
shift at every scattering point. Hence, the phase-angle 
differences of a scattering point within a period are 
distributed over a small angular range and this angular 
range is randomly distributed over [0, 2x]. Since phase 
angle is a cyclic function of 2x, the angular range cannot 
have any 2x jumps when we calculate the phase differ- 
ence of the scattering center of gravity. 


IV. Experimental Results 


In order to confirm the above analysis, a great deal of 
simulation and experimental work was done. In com- 
puter simulation, a flying scattering point model and an 
electromagnetic scattering model were employed. Exper- 
imentally, measured Gata from a rotating model was 
used. Furthermore, simulated flight data was also 
studied. The results are excellent. Nevertheless, we are 
most concerned about imaging of flying targets. 


We were fortunate to obtain a set of Boeing 727 flight 
data from Professor Steinberg of the University of Penn- 
sylvania. Figure 6 shows the image generated after 
motion compensation (aircraft and radar parameters are 
shown in the figure). Comparing to the actual silhouette 
of the aircraft at the bottom, it is not difficult to come to 
the conclusion that the imaging effect is satisfactory. 
This illustrates that motion compensation based on the 
scattering center of gravity is a viable approach. This 
technique will be verified with more experimental data, 
especially actual moving target data, in the near future. 


(a) Sketch of the aircraft 





(b) Aircraft image 
Figure 6. Imaging Based on Motion Compensated Data 





It should be pointed out that the center of mass of the 
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Only in a statistical average sense. Hence, it is necessary 
to smooth out the results while performing complex 
envelope delay compensation. This technique has 
proven to be successful in a variety of experimental 
conditions. As for smoothing of initial phase compensa- 
tion, computer modeling was also successful. However, 
it may not work on actual data because there are other 
factors affecting the initial phase. For instance, insta- 
bility of emission frequency will cause random distur- 
bance of initial phase and this should not be smoothed 
out. 


Motion compensation in ISAR imaging needs further 
research. The method of using isolated points as refer- 
ence points is effective when there are isolated scattering 
points. However, it becomes useless with targets that 
have densely distributed scattering points. In this case, 
the scattering center of gravity method has an edge. 
Furthermore, the more scattering points there are, the 
more even the distribution becomes and the more favor- 
able the compensation ts. In other words, isolated point 
and scattering center of gravity complement each other. 


Acknowledgments 


The authors wish to express their gratitude to Professor 
Steinberg of the University of Pennsylvania for pro- 
viding the actual flight data, to the ISAR research group 
at Harbin Institute of Techn: ogy for providing micro- 
wave darkroom data, and to the 207th Institute of the 
Ministry of Aerospace Industry for supporting our 
microwave darkroom experiment. 


References 


1. D. R. Wehher, “High-Resolution Radar,” Artech 
House, 1987. 


2. Deng Wenbiao and Bao Zheng, “Relation Between 
One-Dimensional and Two-Dimensional Imaging of 
Rotating Targets,” XIAN DIANZI KEJI DAXUE 
XUEBAO [JOURNAL OF XIDIAN UNIVERSITY], 
Vol 17, No 1, 1990, p I. 


3. Deng Wenbiao, “Basic Characteristics of R-D Imaging 
of Rotating Targets,” DIANZI KEJI DAXUE XUEBAO 
(Chengdu) [JOURNAL OF UNIVERSITY OF ELEC- 
TRONIC SCIENCE AND TECHNOLOGY OF 
CHINA], Vol 9, No 1, 1990, p 38. 


4.C. C. Chen and B. C. Anders, “Target Motion-Induced 
Radar Imaging,” IEEE TRANS., Vol AES-16, No 1, 
1980, pp 2-14. 


5. M. T. Pickett and C. C. Chen, “Principle of Inverse 
Synthetic Aperture Radar Imaging,” EASCON, 1980, pp 
340-345. 


Bao Zheng 


Since graduating from Northwest Institute of Telecom- 
municaiions Engineering (presently Xidian University) 
in 1953, he has stayed at the institute to do research and 


aoaaaki, DPacaeamélea, tra ia a aanfeeanoe am tha f ieadssata Gahan 








JPRS-CST-92-023 
19 November 1992 


of Electronic Engineering and an adviser for doctoral 
candidates. His areas of interest include signal pro- 
cessing and radar systems. Professor Bao is a member of 
the Chinese Society of Electronics and an editor of 
DIANZI XUEBAO [ACTA ELECTRONICA SINICA]. 


Deng Wenbiao 


Born in February 1961, he graduated from the Depart- 
ment of Information Engineering at Northwest Institute 
of Telecommunications Engineering in 1982. He earned 
his Master’s in Engineering from the Department of 
Radio Technology at Chengdu University of Electronic 
Science and Technology in 1985 and his Ph.D. in Engi- 
neering from Xidian University in 1988. He is engaged 
in research in areas such as radar imaging, image pro- 
cessing, digital circuitry and digital system diagnosis. 


Yang Jun 


Born in 1965, she graduated from the Department of 
Electronic Engineering of Xidian University in 1987. In 
1990, she earned her Master’s degree in signals, circuits 
and systems from Xidian University. Presently, she is 
working on her Ph.D. degree at the same institution. Her 
thesis is in the area of ISAR imaging and motion 
compensation. 


Motion Compensation in ISAR Imaging for 
Isolated Scatter 


92FE0801C Beijing DIANZI XUEBAO 
[ACTA ELECTRONICA SINICA] in Chinese 
Vol 20 No 6, Jun 92 pp 21-27 


[Article by Meng Jing [1322 7324] and Zhang Zhizhong 
[1728 4160 0022] of Nanjing Institute of Electronic 
Technology, Nanjing 210013: “Motion Compensation in 
ISAR Imaging for Isolated Scatter”; MS received May 
90, revised Feb 91] 


[Excerpts] Abstract 


When there is a relatively isolated and strong scatterer in 
the target, ISAR [inverse synthetic aperture radar] 
motion compensation can be transformed into the align- 
ment of each return and the phase compensation of each 
return. In this work, the return phase is used as an 
approximation to estimate the return Doppler fre- 
quency. In addition, a novel method is presented to 
eliminate any ambiguities of Doppler frequency. The 
target trajectory is re-established using the unambiguous 
Doppler frequency to align the range of each return. 


Key Words: ISAR imaging, motion compensation, Dop- 
pler frequency ambiguity resolution. 


I. Introduction 


In microwave imaging with a phased array antenna, the 
antenna must be aimed at the target or target area for 
focusing. This is to perform phase compensation for each 
element and then find the coherent sum. SAR (synthetic 
aperture radar) or ISAR imaging can be treated as 
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imaging with a large-aperture phased array antenna. 
Therefore, each return should be compensated. The 
difference is that because the aperture is extremely large, 
it is not only necessary to perform phase compensation 
but also range alignment for each return. The combina- 
tion is considered as motion compensation in this paper. 


In studying ISAR imaging, because it is difficult to 
accurately measure the trajectory of the target, motion 
compensation has always been a challenge. In order to 
solve this problem, many assumptions have been postu- 
lated. For instance, the target has an isolated and strong 
scatterer and this small scatterer can be resolved from 
the longitudinal range of the wide transmission signal. 
Furthermore, the scatterer has a small transverse dimen- 
sion (which is less than the required transverse resolu- 
tion). Then, the small scatterer can be treated as a 
relatively isolated source (a point scatterer in an ideal 
situation). It has been demonstrated in reference | by 
means of statistical analysis of a great deal of experi- 
mental data that the probability that such a relatively 
isolated scatterer exists is 80 percent when the range 
resolution is 5, = | m. 


For targets such as airplanes, the isolated scatterers 
might be wing tips, engines and tails. Under the assump- 
tion that such an isolated scatterer could be found, many 
ISAR imaging experiments on real flying targets or 
airplane models have been reported.'~“ 


Motion compensation in ISAR imaging primarily 
involves range alignment for every return. A spatial 
cross-correlation method and frequency-domain method 
were used for range alignment in reference 2. Both 
methods involve a great deal of computation. Here, we 
will introduce a simple and direct method for range 
alignment by estimating the unambiguous Doppler fre- 
quency in order to re-establish the trajectory of the 
moving target. 


Since the radial velocity of the target is not constant, i.e., 
its Doppler frequency is not invariant, the Doppler 
frequency estimated based on the returns accumulated 
over a period At in the past cannot accurately reflect the 
Doppler frequency at the present time. In this work, we 
use the phase of the target return to estimate the Doppler 
frequency of the echo and present a new method to 
eliminate Doppler ambiguities. 


The range alignment method is described below. 


II. Motion Compensation 


In the following discussion, let us assume that the 
operating wavelength of the radar is A = 0.05 m, and the 
range resolution is 5, = | m for the 150 MHz wide-band 
signal emitted by the radar. After taking the computation 
load and the transverse resolution required, the pulse 
repetition frequency (prf) chosen is F, = 100 Hz. 


By using a coherent measurement radar, we can obtain 
data such as the target range R’(t) and complex target 
return a’(t). If discrete sampling is employed, we get: 





12 


Target range R’(i), complex target return a’(ij) (where 1 = 
1, 2, ..., N represents the ith transmission (or reception), 
j= 1, 2, ..., M represents the jth range sampling). 


a’(i,j) can be used to determine the mode phase y’(i,j) of 
the target return phase yw(i,j), where (-2n < y’ < 0). 


When there is a relatively isolated scatterer in the target 
(i.e., an isolated point), it is assumed that the target range 
measured corresponds to the range of the isolated point. 
Assuming that random error in range measurement is 
totally due to thermal noise, then the error €p(i) obeys a 
normal distribution N(O,oRg), where 6g is the rms (root 
mean square) value of the normal distribution. When 
26x = 5p/2, i.€., Gp = 5p/4, the probability of aligning the 
return to the same range unit 5, is 95.45 percent.° 
Because the accuracy of ranging cannot meet the above 
requirement, it is necessary to perform range alignment 
for each return. [passage omitted] 


III. Computer Simulation and Analysis of Results 


With reference to a six-point model, as shown in Figure 
1, computer simulation was done using the method 
described above. The ratio of points (1) through (6) is 
1:1:1:0.5:0.6:0.8. Each point is rigidly connected and the 
points are equally spaced. Let us assume that the target 
height is Y, = 8,000 m and the vertical distance is Zp = 
-39.192 km, then, the nearest range is Ry = 40 km, as 
shown in Figure 2. Assuming that the target is moving 
linearly from point A to B at a constant speed of V, = 250 
m/s, as shown in Figure 3, then AB = V,T = 250 m/s x 4 
s = 1,000 m. The transverse resolution is 


~~ A ~e A ~w 
0s~ Ap DAB/R,  I™ (22) 





Here, the signal-to-noise ratio (SNR) is defined as the 
average SNR of the element of return with the strongest 
average power: 


ifaie: ; 
SNR& =max {D5 1 ai, 4) 1/03} (23) 


ge] 


where | < j < M = 64, N = 400 and a,’ is the standard 
deviation of noise. 


Figures 4-6 show the results obtained at different SNRs. 
The result shown in Figure 6 is considerably worse than 
those shown in Figures 4 and 5. Furthermore, there is 
target movement in the longitudinal direction. This is 
because of an error in the expansion of the ambiguous 
Doppler frequency at lower SNRs which causes the 
linearly varying Doppler frequency to have an abrupt 
discontinuity (up or down), as shown in Figure 7. The 
dotted line in the figure represents the original linear 
Doppler frequency and the solid line indicates the dis- 
continuity. This causes a relatively large error when we 
reconstruct the range and consequently produces more 
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error in range alignment. Of course, this makes it impos- 
sible to produce an accurate image of the moving target. 
As for the cause of the error in Doppler frequency 
expansion, it is due to the fact that the noise-to-signal 
ratio is too high in some instances. There is a large error 
in the ambiguous Doppler frequency derived from the 
phase of the return which leads to an erroneous Doppler 
frequency expansion. Further study on this issue and 
more effective expansion of Doppler frequency will be 
discussed in a later article. 


On the basis of the above discussion, in order to recon- 
struct the trajectory of a moving target, the unambiguous 
Doppler frequency should be accurately estimated first. 
When the target contains a relatively isolated and strong 
scatterer, the key issues in ISAR imaging are how to 
perform range alignment for every return wave and how 
to find the range element for the isolated scatterer. 
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IV. Conclusions 


Under the premise that the target has a relatively iso- 
lated and intense scatterer, this paper presents a way to 
estimate the Doppler frequency of the return by its phase 
and introduces a novel method to resolve any Doppler 
ambiguities. Unambiguous Doppler frequency is then 
used to reconstruct the trajectory of the moving target to 
perform range alignment for every return. Furthermore, 
the isolated scatterer is located based on the least square 
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of normalized return amplitudes. The phase of the 
isolated scatterer is then used as the standard to com- 
pensate for the phase of the return from each range 
element. This method has been found to be effective 
through theoretical analysis and computer simulation. 
However, when such an isolated scatterer does not exist, 
or even when such an isolated point exists but the 
required transverse resolution is high, this phase com- 
pensation method is no longer effective. This problem 
will be discussed in future articles. 
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[Text] Abstract 


It is well known that the range resolution and cross-range 
resolution of a range-Doppler imaging radar are deter- 
mined by the effective bandwidth of the transmitted 
signal and the angle of rotation of the object relative to 
the line-of-sight of the radar during the coherent pro- 
cessing period. The purpose of this work is to lift the 
limitation of conventional FFT [fast Fourier transform] 
range-Doppler processing to significantly improve the 
resolution of a range-Doppler imaging radar. The 
approach adopted is to use linear prediction for data 
extrapolation before performing Fourier reconstruction. 
Preliminary results of rotating platform imaging 
obtained in a microwave anechoic chamber indicate that 
this is a promising technique to enhance the resolution of 
range-Doppler imaging radar. 


Key Words: radar imaging, superresolution processing, 
linear prediction, data extrapolation. 


I. Introduction 


There are a variety of range-Doppler radars for different 
applications. Synthetic aperture radar (SAR) mounted 
on a moving platform is used in air-to-ground imaging. 
Inverse synthetic aperture radar (ISAR) fixed on the 
ground is used in the imaging of airplanes, space objects 
and heavenly bodies. A rotating platform radar is also 
fixed on the ground. It can also be used in imaging 
targets on the rotating platform. These range-Doppler 
radars have a common critical parameter, i.e., resolution 
of image reconstruction. 


A range-Doppler radar gets the required resolution from 
the range-Doppler principle.' The range resolution of the 
image, p,, is dependent upon the effective bandwidth of 
the transmitted signal B: 


p, = C/2B (1) 


where C is the speed of light. On the other hand, the 
cross-range resolution, p,, is dependent upon the angle 
A® through which the object rotates relative to the radar 
line-of-sight (RLOS) in the coherent processing period. 
The relation between p, and A@ is as follows: 


p, = /2A0 (2) 
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where A is the radar wavelength. It should be pointed out 
that equations (1) and (2) are rigorously valid only when 
the A®@ is so small that the scattering point of the target 
does not go over a range resolution unit and a Doppler 
resolution unit. With classic signal processing tech- 
niques, it is impossible to overcome this limitation 
imposed by equations (1) and (2). In this case, the 
effective bandwidth of the transmitted signal and the 
total angle of rotation of the target with respect to the 
RLOS must be increased to enhance the range resolution 
and cross-range resolution of the image, respectively. It 
is often not desirable to do so in many applications. 
Another approach to improving the resolution capability 
is to use advanced, rather than conventional, signal 
processing methods because these methods may be able 
to enhance the range resolution to a level that far exceeds 
the limitation discussed earlier. 


Signal processing in range-Doppler imaging can be con- 
sidered as a two-dimensional spectral analysis problem. 
Many high-resolution spectral estimation and frequency 
estimation methods, i.e., modern spectral estimation 
methods, have been studied.” In order to take advantage 
of the superiority of these methods, an accurate signal 
model and a sufficient signal-to-noise ratio (SNR) are 
required. Furthermore, modern spectral estimation 
methods are far more complex compared to the classical 
Fourier method. For this reason, people continue to 
investigate methods that combine modern high- 
resolution techniques with the classical approach. One of 
the hybrid approaches? is relatively simple in concept 
and in computation. In this approach, a linear prediction 
filtering method is employed to extrapolate the data 
beyond the observation window. Then, the extrapolated 
data undergoes discrete Fourier transform (DFT) to 
attain a higher frequency resolution. This method is used 
to improve the resolution of rotating platform imaging to 
construct a two-dimensional image of a B-52 model by 
processing data obtained in a microwave anechoic 
chamber. Preliminary results indicate that this superres- 
olution range-Doppler radar imaging method is less 
complex in computation and more general in applica- 
bility. In the area of modern spectral estimation methods 
to improve radar imaging resolution capability, 
Gabriel*> introduced a maximum likelihood filter 
method and a linear prediction filter method to achieve 
superresolution range-Doppler imaging. These methods 
involve the determination of the inverse of matrices, 
which requires a great deal of computation. Further- 
more, references 4 and 5 only provided imaging results 
of simulated data and did not verify the methods with 
measured data. 


In this paper, the use of linear prediction for data 
extrapolation to achieve superresolution imaging 1s 
described first. A similar method is then used to obtain 
superresolution of cross-range imaging. Finally, a super- 
resolution two-dimensional image of a B-52 model is 
presented. 
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II. Superresolution Range Imaging 


The high range-resolution capability needed in range- 
Doppler radar imaging can be realized by using tech- 
niques such as bandwidth pulse compression, time- 
bandwidth interchange, or frequency-stepping 
waveform. The method described here is also applicable 
to pulse compression and time-bandwidth interchange. 
The difference is that the pre-treatment is slightly 
different.'-° 


When frequency-stepping waveform is employed, the 
radar transmission frequency and the receiver local- 
oscillation frequency vary in a stepwise manner. Figure | 
shows a target at a distance R away from the radar; f(r) is 
the range function of the reflectivity of the target. It is 
also known as the range image of the target. The variable 
r represents the distance of the target with respect to the 
reference point o. With a specific viewing angle, the 
radar transmits a series of K pulses at a frequency 


f, = f, + kéf, k = 0, l, sovcce ,K- 1. (3) 


The return from the target after compensating for the 








Sir) 
Radar Target:/~ 
R——_—_ 








phase term due to range R can be expressed as follows: 





$.= AU P07 exp {+i & farhdr, 


(4) 


where A is a proportionality factor which is dependent 
upon factors such as power and signal transmission. It is 
also assumed the pulse width covers the entire target 
range. Integration is carried out along r over the entire 
target range. Equation (4) implies that S, is the inverse 
Fourier transform of f(r). Hence, after knowing the K 
vector of S,, it is possible to estimate f(r) by computing 
the DFT of S, where k = 0, 1, ...., K - 1. Obviously, the 
range resolution of such a target image is 


p, = C/2K8f (5) 


This is consistent with equation (1) because here the 
effective bandwidth of the transmitted signal is B = K8f. 


Let us assume that the scattering from an extended target 
comes from a finite number of scattering centers.’ For an 
individual scattering center, if its range is fo, the 
observed signal in the frequency domain is a complex 
sine function of 











Af (19) exp{ + j x ref} 


C 


Therefore, the observed signal of the target is comprised 
of a series of complex sine functions in the frequency 
domain. It can be accurately predicted by a linear 
prediction filter of an order equal to or greater than the 
number of complex sine functions.” Although this state- 
ment is valid only when there is no noise in the observed 
signal, in reality it is still possible to employ linear 
prediction to reliably extrapolate the data outside the 
observation window when the SNR is relatively low. 


Assuming that actual data S,, k = 0, 1, ...., K- 1, has been 
obtained at a given viewing angle, a single-step linear 
prediction filter of order P is used to extrapolate the 
data. The filter coefficient is a,, p = 1, ...., P. This filter 
produces a forward extrapolation sample S, as follows: 


S,= DS 08 k>K (6) 


P=1 


The backward extrapolation process is similar. The 
difference is that the filtering direction is reversed. 
Moreover, the filter coefficients are complex conjugates 
of those of the forward filter. In other words, the extrap- 
olation takes place according to the following: 


P 
Sp =D S 44993, K<1 (7) 
?=1 


Thus, the data is actually measured if | < k < K and the 
data is extrapolated when k > K and k < |. The forward 
extrapolation length is equal to that of backward extrap- 
olation. The Burg maximum entropy (ME) technique is 
used to calculate coefficients a,, p = I, ...., P. A vector 
with N data points (N > K) can be obtained with the 
method described above. This N-point extrapolated 
vector replaces the original K-point vector in the DFT 
computation. The range image thus derived can be 
expected to have a much better range resolution. 


Figure 2 shows the overall range image of a B-52 model. 
This is the result of computer processing based on 
microwave anechoic chamber data. A frequency step- 
ping technique was used in the experiment. The number 
of frequency steps is K = 32. The transmitted signal has 
an effective bandwidth of B = 2.15 GHz. Table | shows 
the raw data obtained in the microwave anechoic 
chamber at a specific viewing angle. Figures 2(a) and (b) 
are the range images obtained with the original 32 data 
points and with 128 extrapolated points with a Ham- 
ming window, respectively. It is apparent that data 
extrapolation by linear prediction, in conjunction with 
conventional Fourier transformation, can provide super- 
resolution range images. 
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Figure 2. Range Image of B-52 Model at «Specific Viewing Angle (a) with 32 original points, (b) with 128 extrape- 
ted points 








Table 1. Raw Data at a Specific Viewing Angle 


















































-.080 99 “11520 -.12826 10653 17118 07572 ~.19275 00786 
-.169 40 -.08068 -.13269 ~.18262 -,06945 ~.26023 -.01409 27528 
048 96 21143 08769 -.08913 09758 06676 08207 19837 
026 57 26618 -.04690 25384 ~.11979 18449 17199 06819 
-.184 57 -.05550 ~.16452 14741 ~.10523 19206 ~.02174 ~.19634 
064 16 ~17319 13968 ~.12860 19338 08315 21725 -.03815 
21441 01149 16985 06559 11442 13303 03429 17366 
015 34 17477 -.07138 13274 -.11886 03660 ~.14958 -.06196 
III. Superresolution Cross-Range Imaging On the basis of the same argument to attain super range 


In small-rotation-angle imaging with the FFT range- 
Doppler process, range and cross range are processed 
separately. As described in Section II, we can obtain a 
range image at every viewing angle. At M equally spaced 
viewing angles, M range images of the target can be 
obtained. Thus, a data vector of M points can be 
obtained with respect to each range unit. The cross-range 
image of a given range unit, i.e., the reflectivity of the 
target in that range unit as a function of the cross range, 
can be derived from the DFT of this M-point vector. 
When it is applied to a sample of M actually measured 
data points, the cross-range resolution is 


p, = A/2M80 (8) 


where 80 is the spacing of two neighboring viewing 
angles. Since A@ = M80, equation (8) is consistent with 
equation (2). 


resolution discussed in Section II, superresolution cross- 
range images can be obtained by using linear prediction 
for data extrapolation in conjunction with conventional 
Fourier transformation. Figure 3 shows the overall cross- 
range image of the B-52 model at a range unit derived 
from measured data. The total angle of rotation of the 
airplane model in the imaging process is A® = 15° and the 
radar wavelength is A = 2 cm. Table 2 lists the 32 original 
data points corresponding to the range unit. Figures Xa) 
and (b) show the result obtained from the original 32 
data points and that derived from the 128 extrapolated 
points with a Hamming window, respectively. It is 
apparent that using linear prediction for data extrapola- 
tion can also significantly improve the cross-range reso- 
lution of range-Doppler radar imaging. 
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Figure 3. Cross-Range Image of B-52 Model. (a) with 32 original points, (b) with 128 extrapolated points 
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Table 2. Raw Data at a Given Range Unit 





















































“391 29° 43871 -.07534 -.23074 50295 -.80743 53297 65919 
-,569 85 -.42205 -1.20443 -.69226 -.40890 -.36944 21238 43367 
303 96 82409 19226 46031 17158 -.13421 11666 61985 

. 097 28 -.10557 -.27174 52000 -.16635 -.05365 10964 39470 
395 75 54619 62743 19806 $2085 -.28259 75538 51553 
319 49 -.41673 -.17068 -.27398 -.13106 -.10458 21626 16420 
479.72 42178 43015 29670 23022 -.15925 01776 -.47091 
-.237 82 -.43817 -.51430 -.18886 -.53399 -.04538 -.26843 02514 

_ IV. Superresolution Two-Dimensional Imaging Figure 4(b). 


The range resolution and cross-range resolution in two- 
dimensional range-Doppler imaging can both be improved 
with the method described above. The results are shown in 
- Figure 4. Using the data shown in two previous sections, a 
two-dimensional radar image of the airplane model was 
obtained by processing the range dimension first for every 
viewing angle and then processing the cross range for every 
range unit. Figures 4(a) and (b) show the images of the B-52 
model processéd based on the original 32 x 32 data matrix 
and the extrapolated 128 x 128 matrix, respectively. A 
comparison shows that our approach significantly improved 
range resolution and cross-range resolution of the image. 
Compared to the photograph of the B-52 model shown in 
Figure 5, the two engines under the wing were clearly visible 
‘in the superresolution two-dimensional image shown in 





Figure 4. Two-Dimensional Image of B-52 Model. (a) 
32 x 32 true data matrix, (b) 128 x 128 extrapolated 
data matrix 








The microwave anechoic chamber data used has a very high 
SNR. It almost can be treated as a noise-free case. In order 
to investigate the effect of SNR on the superresolution 
imaging technique described, white Gaussian noise was 
added to the microwave anechoic chamber data. A great 
deal of results indicate that this method can still produce 
better superresolution two-dimensional images even with a 
low SNR. Figure 6 shows the results with an SNR of 0 GB. 





Figure 6. Two-Dimensional Images of B-52 Model With 
SNR = 0 dB. (a) 32 x 32 true data matrix, (b) 128 x 
128 extrapolated data matrix 





V. Conclusions 


Our preliminary results indicate that superresolution 
range-Doppler radar imaging using data extrapolation 
via linear prediction is a relatively simple and practical 
superresolution imaging method. It is capable of lifting 
the resolution limitation of conventional Fourier image 
processing. This method can produce higher resolution 
images with identical effective transmission bandwidth 
and total angle of rotatio.1 during imaging. Or, it can 
produce an image of sim lar resolution with a smaller 
bandwidth and total rotation angle. 


During the process of revising this paper, the authors 
found that Li® et al. used a similar method to improve 
the resolution in one-dimensional range imaging by 
extrapolating data in different isolated segments. This 
work has presented a method for two-dimensional super- 
resolution range-Doppler imaging. In addition, the effect 
of SNR on this method has also been studied. Even with 
a low SNR, this method was still found to be capable of 
producing better superresolution two-dimensional 
images. 


The authors wish to express their gratitude to Professors 
Liu Yongtan [0491 3057 0962], Cao Zhidao (2580 1807 
6670], and Meng Xiande [1322 2009 1795] of Harbin 
Institute of Technology and to colleagues at the 207th 
Institute of the Ministry of Aerospace Industry for 
providing the microwave anechoic chamber data. 
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[Text] Abstract 


The imaging of targets rotating in a small azimuth angle 
range is discussed. A unified image reconstruction 
scheme suitable for bistatic, monostatic, near-field and 
far-field imaging geometries is derived. Super-resolution 
imaging techniques are investigated. A new imaging 
reconstruction algorithm, combining spectral extrapola- 
tion and coarse “shape’’ information, is presented. 
Finally, experimental results in microwave anechoic 
chambers are described. 


Key Words: microwave imaging, radar targets, imaging 
processing. 


I. Introduction 


Microwave imaging of rotating targets has been widely 
used in many areas, such as: (1) Target diagnosis and 
identification. Rotating target imaging is a basic tool for 
studying the electromagnetic scattering of targets. It 
plays an important role in the study of electromagnetic 
scattering mechanism and RCS reduction for stealth 
research and target identification. (2) SAR and ISAR 
signal processing. SAR or ISAR imaging with motion 
compensation is equivalent to rotating target imaging. 
(3) Non-destructive testing. The principle of rotating 
target imaging is also widely used in non-invasive and 
non-destructive testing. Therefore, in-depth study of 
rotating target imaging and related processing tech- 
niques is a basic research project with a wide range of 
applications. 


D. L. Mensa! once pointed out that imaging of targets 
rotating over 360° could achieve extremely high nominal 
two-dimensional resolution (a fraction of a wavelength 
of the radar). The condition to attain such a high 





resolution is that all the scatterers on the target must be 
isotropic. Furthermore, the scatterers must not move in 
order to obtain a uniform synthetic radar aperture over 
360°. However, the following two reasons make rotation 
imaging in small azimuth angle range more practical. (1) 
Real radar targets often cannot satisfy the requirement 
that the scatterers be isotropic and stationary. Hence, it 
is not possible to attain the high theoretical nominal 
resolution. Furthermore, the amount of work involved in 
data acquisition and processing for 360° rotating target 
imaging is tremendous. (2) As radar signal technology 
advances, high radial resolution of the target can be 
obtained with a wideband signal alone. It usually takes a 
small azimuth angle to achieve high transverse resolu- 
tion. Therefore, imaging of targets in a small azimuth 
angle range has been developed for a variety of applica- 
tions in the past two decades. 


II. Image Reconstruction Algorithm 


The bistatic geometry for rigid target imaging is shown in 
Figure |. The target coordinate is (x-y) and the radar 
coordinate is (u-v). 


x = u cos 6 - v sin 8, y = u sin 6 + v cos 
8 (1) 


or 


x sin 8 + y cos 


u = x cos 6 + y sin 8, v=- 
8 (2) 











¥% d =~ 
Transmission Receiving 
antenna antenna 


Figure 1. Geometry for Microwave Imaging 








The two-way range between a poirt on the target (x,y) 
and the transmitter and receiver is 


R,=CCR,+¥)?+ (d+u)?}? 
= {Ri+7?7+4'+2r7R, sin (9—§) 
+dcos (@—§)]}” 
R,=((R,+v)?+ (d—u)?)” 
= {R2+r°+47+2r[R, sin (9 —§) 
—dcos (9 —9)]}*? 


(3) 


(4) 


and 
R=R,+R, 
= {Ri+r+a+2r[R, sin (9—§) 
+dcos (9—§)}}!? 


+{R3+r+24+2r[R, sin (9—§) 


. —dcos (P—9) J} (5) 
From the above formulas, after frequency mixing, cali- 
bration and filtering, the return can be derived as fol- 
lows: 


F(K,6)=|["_ f(r, #)exp{—74K(R, +R, 


—2R,)} rdrd¢ (6) 
From equation (6), the return is the correlation integral 
of the scattering signals along various equidistance 
lines; f(r,@) is the scattering function of the target. In 
two-dimensional imaging of rotating targets, the spec- 
tral domain signal of f(r,@) is measured in the ring- 
shaped spectral domain shown in Figure 2. The inte- 
gration for image reconstruction is also performed in 
this domain. The image reconstruction equation 
expressed in polar coordinate is as follows: 


Jerson [Tf FRO Rexp (iA KCR, +R. 
—2R,)} 4X48 (7) 








Figure 2. Polar Coordinate Data Format in Imaging of 
Rotating Target 
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Various reconstruction methods can be designed based 
on equation (7), including the two-dimensional Fast 
Fourier Transform (2-D FFT), filtered back-projection 
(FBP) and Fourier domain extrapolation. 


Experimental and theoretical studies indicate that the 
FBP method not only is more flexible but also produces 
higher quality reconstructed images in two-dimensional 
imaging of rotating targets. It is suitable for anechoic 
chamber imaging. This method is described as follows: 


In equation (7), if the integration of K is treated as the 
inverse Fourier transform of variable |, = (R, + R, - 
R,)/2, then the reconstructed f[caret}(r,o) becomes 


K2-Ki 
P.)= |) IK+K,|F(K+Ky,6) 


a oy? 
fer, 9)=["F Po Cl sdexp(i27Kyl )¢6 


1-=(R,+R,—R,)/2 (8) 


In actual computation, the first formula (radial focusing) 
is computed by FFT and then the second formula is used 
to compute back-projection (transverse focusing). It is 
very similar to the FBP technique in tomographic 
imaging. Hence, it is also called filtered back- 
projection.” 


The third formula in equation (8) describes the path of 
back-projection. From analytical geometry, the path of 
back-projection described by that formula follows an 
elliptical trajectory. It is identical to the integration path 
of radar projection measurement, i.e., 


1,={[Ri+°+2?+2r(R,sin (P—§) 
+dcos (9—Q)) }? 
+(R5+1r7+2?+2r(R, sin (9—§) 


—dcos (9—9)) }"—2R,} /2 (9) 


If it is a Monostatic measurement, then in Figure | d = 0 
and equation (9) becomes 


l, = [R? + r? + 2r Ro sin (@ - 8)]'”? - Ro (10) 


This describes a family of circles. In other words, in 
monostatic measurements, the projection integration 
path is a circular arc and the back-projection path in 
image reconstruction is also circular. 
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Furthermore, when R, >> r, i.e., the far-field condition 
is met, then 


l. = r sin (@ - 8) (11) 


Equation (11) describes a straight line. It is suitable in 
far-field imaging (plane wave front). 


The derivation of equation (8) does not require the radar 
to emit plane electromagnetic wave fronts. Therefore, 
the resultant reconstruction method is applicable to a 
varicty of imaging conditions such as monostatic, 
bistatic, far-field and near-field. Incidentally, because 
the target spectrum obtained in microwave imaging of a 
rotating target is located in a ring-shaped domain, the 
expansion function of the imaging system has very high 
side lobes. Therefore, usually a smoothing window (at 
small azimuth angle) or variable trajectory filter (360° 
imaging)° is required in the actual reconstruction com- 
putation to suppress the side-lobe voltage in the recon- 
structed image to improve image quality. 


III. Super-Resolution Imaging Technique 


It is well known that linear prediction is widely used in 
the estimation of high-resolution power spectrum of a 
discrete time sequence, especially when the sequence is 
very short. By means of extrapolating the sequence, it is 
possible to achieve higher resolution than that of a 
conventional Fourier spectrum. Since in microwave 
imaging the two-dimensional resolution is determined 
by the frequency spectrum window and azimuth 
window, when the radar bandwidth is narrow or the 
azimuth angle is small, the resolution deteriorates rap- 
idly. In order to obtain a high-resolution microwave 
image with a small amount of data, linear prediction is 
used to investigate the super-resolution imaging of 
microwave data. 


Based on the Burg algorithm,®’ a sequence with M 


specimens [X,, Xz, X3, .--, Xm] iS extrapolated by linear 
prediction to have N (N > M) specimens. The predic- 
tion parameters can be obtained by minimizing the sum 
of errors in energy, E,, of forward and backward predic- 
tion, 1.€., 


N-1 N-1 
Eo= Dil Eo: l'+ Sil Evol (12) 
« =?7 


where E,,, is the forward prediction error of the pth order 
model: 


a 
Eny=> jo onre-e 
ke 


(13) 


and E 
model: 


is the backward prediction error of the pth order 


(14) 


? 
—_—_ - 
Eos = D523 eX e-9-k 


a=d 








Here, a,, is the linear prediction parameter and * repre- 
sents a complex conjugate operation. 


The linear prediction equation for extrapolation with M 
specimens is as follows: 


bid 
us j = D> eetu-j-k 7>0 (18) 
k=l 
ad 
BO X58 1>0, (16) 
k= 


where [caret] represents the estimated value. 


Figure 3 shows schematic diagrams of one- and two- 
dimensional extrapolation by linear prediction with par- 
tial data in the target spectrum. After spectral expansion, 
the image can be reconstructed with the resultant data 
using the method described in the last section to obtain 
a super-resolution image of the target. 


Azimuth Frequency 
spectral spectral 
expansion expansion 


(a) One-dimensional spectral 
expansion 


Measured 
CO Predicted 


(b) Two-dimensional spectral 
expansion 
Figure 3. Image Spectrum Expansion by Linear 
Prediction 





In image reconstruction with extrapolated data by linear 
prediction, it was found that although the resolution is 
enhanced, the side-lobe voltage of the image is also 
increased due to prediction errors. In order to suppress 
image side lobes, an algorithm is developed to combine 
linear prediction with the rough “shape” of the target. 
The “shape” is derived from the original data. Specifi- 
cally, both original data and extrapolated data are pro- 
cessed in image reconstruction. In radial focusing with 
FBP, two sets of range distribution plots are obtained. 
The range distribution plot obtained from the extrapo- 
lated data is compared to that derived from the original 
data to eliminate the side lobes. After both images are 
reconstructed, comparison can be made again to further 
eliminate any side lobes other than the “shape”’ of the 
target in the original image. This new algorithm can 
drastically lower the side-lobe voltage of the image. 
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Figure 4 shows the flowchart of this method. The results 
will be presented in the next section. 


| “ghase* data ~4 


|Frequency spectrum expansion! 


spectrum expansion 
1 


[Image recone truction} 


Target image!) 
Figure 4. Flowchart of the New Image Rev onstruction 
Algorithm 
































IV. Anechoic Chamber Experiments and Results 


The entire experiment was done on the microwave 
measurement and imaging system of an HP8510 net- 
work analyzer. The major specifications of the system 
are shown in Table |. 


Table 1. Major Specifications of the Microwave Imaging 


System 
Sweep range 2-4, 4-8, 8-18 GHz 
Frequency stability 1x 10°9/day 
Output power 20 dBm 
Antenna gain 23 dB 
Rotating platform accuracy 0.05° 
Minimum step of rotating plat- 0.1° 
form 
Amplitude measurement accu- +/- 0.5 dB 


racy 


Sweep linearity 1% in sweeping mode 
Sweep mode Stepping, sweeping 
Receiver sensitivity -85 dBm (unaveraged) 
-103 dBm (averaged over 64 
times) 
Dynamic range 75 dB (unaveraged) 
93 dB (averaged over 64 times) 
Phase measurement accuracy +/- 3° 


During the imaging process, a microcomputer controls 
the transmitter, receiver and azimuth angle. The ampli- 
tudes and phases of the wide-band scattered returns from 
the target are collected at different azimuth angles. The 
data gathered is also pre-processed by the microcom- 
puter. 


The pre-processing, including software range gate tech- 
nology for background cancellation, data scaling and 
clutter filtering outside the target range, can be expressed 
by the following equation: 


ae Daye (17) 


C=E-Fts—py PS 
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where E is the accurate amplitude phase of the cali- 
brator, T is the frequency-domain measurement of the 
target, T, is the intrinsic scattering of the target, S is the 
frequency-domain measurement of the calibrator, S, is 
the inherent scattering of the calibrator, B, is the back- 
ground during target measurement, B, is the background 
during scalar measurement, and script F[ ] represents the 
filtering process. 


The RCS amplitude of the pre-processed imaging data 
can be expressed in dB referenced to a square meter 
(dBsm). The phase is centered around the rotating plat- 
form. The data is then used for image reconstruction and 
image Output by a microcomputer. 


Figure 5 shows the resultant image of a reduced B52 
model in an anechoic chamber. The experimental con- 
ditions are: center frequency | 1.5 GHz, radar bandwidth 
7.0 GHz, frequency step 35 MHz, target rotation angle 
20°, angular sampling interval 0.5°, antenna polarization 
state H-H, test range 12.5 m. A Hamming window was 
added in image reconstruction. The specifications 
attained are: radial resolution approximately 4 cm, 
transverse resolution 7 cm, and image side-lobe voltage 
less than -40 dB. The dynamic range of the imaging 
system is actually limited by the side-lobe voltage. The 
dotted line in the figure represents the actual silhouette 
of the reduced model. Its scattering center and the 
contour lines obtained from imaging are in excellent 
agreement. 


Figures 6 and 7 show the reconstructed images derived 
from the raw data and the two-dimensional extrapolated 
data, respectively. Some data were acquired with a 
reduced bandwidth of 3.5 GHz and smaller azimuth 
angle of 15°. Comparing Figure 7 to Figure 5, we find 
that the resolution of the reconstructed image from the 
raw data and that of the reconstructed image based on 
extrapolated data are essentially the same. However, 
there is some error in the amplitude of the scattering 
center (approximately 1-2 dB). 


V. Conclusions 


(1) Using a 7 GHz wide-band frequency spectrum and 
20° rotation angle aperture imaging, it is possible to 
obtain a high-resolution microwave image of the order of 
a centimeter. 


‘2) The FBP algorithm derived is applicable to a variety 
of imaging geometries and conditions including bistatic, 
monostatic, near-field and far-field. 


(3) By combining spectral expansion and target rough- 
shape data, super-resolution imaging can be attained. 
Consequently, the amount of imaging data required is 
drastically reduced. At 3.5 GHz and 15° azimuth 
window, the resolution is almost as high as at 7 GHz and 
20° azimuth window, and side-lobe voltage is better than 
-40 dB. However, there is some error in amplitude. 
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Figure 5. Imaging of a BS2 Model 
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Figure 7. Image Reconstructed From Spectral Expanded Data 
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Finally, in addition to spectral expansion, the quality of 
the image may be further improved by using polarization 
information and target symmetry. This effort is in 


progress. 
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[Excerpts] Abstract 


In this work, the scattering of a large planar slot-array 
antenna is studied using the momentum method, phys- 
ical optics method and singularity elimination tech- 
nique. The radar cross section (RCS) of a typical array 
was calculated and verified experimentally. It is found 
that the calculated result is in good agreement with the 
experimental data. The computations also show that the 
contribution of the plane to the scattering field (structure 
term) is far greater than that of the slot (mode term). 


Key Words: planar slot-array antenna, radar cross sec- 
tion, scattering. 


I. Introduction 


The large planar slot-array antenna is one of the com- 
monly used airborne antennas. However, it produces a 
very strong RCS contribution in the direction of the nose 
of the aircraft. In order to improve the stealth character- 
istics of the aircraft, it is necessary to understand the 
scattering mechanism of this type of antenna and ana- 
lyze and calculate its RCS contribution. To date, most of 
the studies on the coupling of incident wave and slot- 
array antenna, or scattering computation, involved the 
combination of a single slot and a cavity.'~* There is no 
information available on the scattering mechanism of a 
large planar slot-array antenna and the detailed analysis 
of its RCS. This work provides a theoretical analysis of 
this type of antenna by rigorously taking the interaction 
between elements into consideration using momentum 
method, physical optics method and singularity elimina- 
tion technique. Moreover, the computing time is drasti- 
cally reduced through the use of symmetry technique, 
Toeplitz property and an excellent inverse technique. 


II. Antenna Array Analysis 


Figures | and 2 show the antenna array. Let us assume 
that the incident direction is on the yoz plane. The 
number of slot vibrators along the x-axis over a length of 
21 is M, the distance between slots is Ax, the slot width is 
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W, the number of vibrators long the y-axis is N and the 
spacing is Ay. The contribution of the entire antenna 
array to the RCS is the phase sum of two terms, i.e., the 
mode term and structure term.‘ 








Array Antenna 


Figure l. 


Zone a: Zone b 


onductor 
Incident plate 
direction 
al y 
Figure 2. Side View of Array 





1. Computation of Mode Term 


The plate can be considered infinite. Its conductive wall 
divides the space into zones a and b. Using the equiva- 
lence theorem [See Reference 5, Sec. 3-5], the problem 
can be split into two, as shown in Figure 3. When the 
aperture is blocked by the conductive plate and replaced 
with an equivalent magnetic flux M, [I = (m-1)N + n, m 
= 1,2,...,M,n=1, 2, ..., NJ, the field in zone a remains 
unchanged. When the aperture is blocked by the conduc- 
tive plate and replaced by an equivalent magnetic flux 
-M,, the field in zone b remains constant as well. [passage 
omitted] 


III. Results and Analysis 


As an example, computation was done on a 4 x 8 
half-wavelength slot-array antenna with A = 3.2 cm. The 
slot width is W = 0.0125A. The end-point spacing in the 
x direction is a half-wavelength and the slot spacing in 
the y direction is also a half-wavelength. The plate size is 
4.06254 x 10A. Figure 4 shows the RCS of the antenna as 
a function of observation angle 8. The unit of RCS in the 
figure is dBA”. From the figure one can see that the 
mode-scattering term generated by the excitation of the 
incident wave at the slot is far less than the structure 
term when the mutual coupling between slots is taken 
into consideration. For instance, it is 21.340 dB less at 
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Figure 3. Equivalent Zones 





the axial direction. Hence, the resultant f° ~ scattering 
field essentially varies according to the patiern of the 
plate structure term in the range of 0-60°. Therefore, it is 
quantitatively demonstrated that structural scattering 
from the plate is the dominating term in the mechanism 
of a slot-array antenna. In order to effectively reduce the 
RCS of a slot-array antenna, it is necessary to minimize 
the contribution of the plate. For example, a traveling 
wave slot-array antenna (non-resonant) may be used to 
offset the peak of its emission from its peak direction of 
scattering to lower the RCS at the nose of the aircraft. In 
addition, a conformal curved surface, such as a cylin- 
drical, conical, spherical, ellipsoidal or paraboloidal sur- 
face, can be used instead of a flat plate to significantly 
reduce the contribution of the structure term to the 
RCS. 





-—-- Plate term 


———Slot term 
Total scattering 


Figure 4. RCS of Planar Slot-Array Antenna vs. Obser- 
vation Angle 





An experimental antenna was made based on the param- 
eters of the model used in the computation. The RCS 
was measured with an outdoor quasi-monostatic RCS 
measurement system (see Figure 5). The measured result 
is compared to the theoretical value, and they are found 
to be in good agreement. The rotation angle a corre- 
sponding to the maximum RCS value in Figure 5 should 
correspond to 6 = 0° in Figure 4. 
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Figure 5. Measured RCS of Planar Slot-Array Antenna 





IV. Conclusions 


Theoretical analysis and numerical computation of the 
RCS of a planar slot-array antenna was carried out using 
the momentum method, physical optics method and 
singularity elimination technique. It was found that the 
structure term plays a dominant role in the scattering 
mechanism. Viable technical approaches of engineering 
value are presented to reduce RCS. 
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[Excerpts] Abstract 


The backscattering of an electromagnetic missile with its 
energy decaying as r* by an ideal metallic elliptical 
cylinder and a sphere is described. The study shows that 
the rates of energy decay for the backscattered fields 
from the elliptical cylinder and sphere are proportional 
to r¢'*® and r¢?*®), respectively. 


Key Words: electromagnetic missile, backscattering, 
backscattered field. 


I. Introduction 


The concept of an electromagnetic missile was first 
presented by T. T. Wu in 1985. It is a transient electro- 
magnetic wave of finite energy that does not decay 
according to r? as a conventional spherical wave in the 
far field. Instead, it decays slowly at an arbitrary rate. 
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This excellent slow decay characteristic of the electromag- 
netic missile can be used to increase the effective range and 
improve the performance of radars.' Properties such as 
short duration (picoseconds) and high peak power (giga- 
watts) can also be used to enhance the detection capability 
of the radar for low-RCS or stealth targets. Therefore, there 
is a need to have a profound understanding of the genera- 
tion and propagation of electromagnetic missiles from both 
theoretical and practical viewpoints. 


In the past few years, a great deal of theoretical and 
experimental work has been done on electromagnetic 
missiles. Several electromagnetic missile launchers were 
constructed. This demonstrated that it is not that diffi- 
cult to make electromagnetic missiles a reality. Besides 
constructing a variety of electromagnetic missile 
launchers, there is another research topic of considerable 
significance, i.e., the study of scattering characteristics of 
electromagnetic missiles from different target shapes. 


It is well known that the energy of a conventional electro- 
magnetic wave decays according to r'' in two-dimensional 
space. Hence, the energy of scattered field decays as a 
function of r°. In three-dimensional space, the energy of a 
conventional electromagnetic field decays as a function of 
r? and that of the corresponding scattered field decays 
according to r*. Recently, Myers and T. T. Wu have 
studied the backscattering of electromagnetic missiles 
from a plate. The results show that the energy of the 
backscattered field of an electromagnetic missile decaying 
as a function of r* from a flat plate still decays according 
to r®(0 < € < 1). This is much slower than the decay of 
the energy of a conventional scattered field. 


This study involves the backscattering of electromag- 
netic missiles from ideal elliptical cylindrical and spher- 
ical conductors of infinite length. The results indicate 
that the backscattered fields of the elliptical cylinder and 
sphere decay according to rt'*® r‘?*®, respectively. 
They are considerably faster than the decay of the 
backscattered field from a flat plate. However, they are 
slower than the decays of conventional two- and three- 
dimensional scattered fields, which are proportional to 
r? and r“, respectively. [passage omitted] 


IV. Conclusions 


This paper describes the backscattering of an electro- 
magnetic missile decaying as a function of r* from an 
elliptical cylinder conductor and a sphere conductor. An 
analysis shows that their scattered fields are proportional 
to r¢'*®) and r(?*®), respectively. This conclusion is an 
important reference point in counter-stealth technology. 
As far as counter-stealth is concerned, an electromag- 
netic missile is an excellent electromagnetic pulse. When 
used in impulse radars, it can effectively counter the 
effects of [radar] absorbent materials (RAM) and profile 
stealth technology. The signal-to-noise ratio of the radar 
can be improved and the effectiveness of the radar can 
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be enhanced. It is of great significance to further inves- 
tigate the theory and practice of backscattering of elec- 
tromagnetic missiles from other targets. There is a great 
deal of work still to be done in this field. 
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